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Abstract 

We calculate the ratios of B and D meson decay constants by applying the 
variational method to the relativistic hamiltonian of the heavy meson. We adopt 
the Gaussian and hydrogen-type trial wave functions, and use six different poten¬ 
tials of the potential model. We obtain reliable results for the ratios, which are 
similar for different trial wave functions and different potentials. The obtained 
ratios show the deviation from the nonrelativistic scaling law, and they are in a 
pretty good agreement with the results of the Lattice calculations. 


TAGS number(s): 12.39.Ki, 12.39.Pn, 13.20.Fc, 13.20.He, 13.25.Ft, 13.25.Hw, 
14.40.Nd 


1: e-mail: dshwang@phy.sejong.ac.kr 
2: e-mail; gkim@cmp.sejong.ac.kr 


The knowledge of the decay constant of the B meson fs is very important, since 
it affects the magnitude oi B — B mixing and the size of CP violation significantly. 
There have been intensive theoretical and experimental researches for improving 
its understanding. However, its theoretical calculation is difficult because it is in 
the realm of nonperturbative QCD and the motion of the light quark in B meson 
is relativistic. Understanding the decay constant better is also invaluable because 
its information can reveal the inside structure of the hadron. 

Grinstein [Q observed that the double ratio of the decay constants {/bs/Ib^)/ 
ifos/foa) is very close to 1 with small correction of order ms/rriQ. He calculated 
the double ratio with the heavy quark effective theory, and obtained 0.967. He 
remarked that the value of which is an important factor for the relative 

strengths of Bg — Bg and Bd — Bd mixings, can then be given reliably from the 
knowledge of the measurable fos/fn^- Oakes [Q also calculated the double ratio 
based on the assumption that chiral symmetry is broken by quark mass terms in 
the Lagrangian. He obtained its value as 1.004, and emphasized the importance 
of the fact that this double ratio is very close to 1. 

When one treats the heavy-light meson in analogy with the nonrelativistic sit¬ 
uation, one expects the scaling law /b//d — \JMd/Mb, since the reduced masses 
of the light quark {u or d quark) in B and D mesons have almost the same value, 
and fpMp = 12 1-0(0) p by the Van Royen-Weisskopf formula @ for the pseu¬ 
doscalar meson P, where -0(0) is the wave function at origin of the relative motion 
of quarks @]. However, the light quark inside B or D meson has large velocity, 
and its nonrelativistic treatment is not legitimate. Indeed, our calculation of the 
decay constants in the relativistic quark model, which we present in this Letter, 
shows that the nonrelativistic consideration is much deviated by the relativistic 
motion of the light quark, since this relativistic nature makes the 0(0) of B and 
D mesons different appreciably. This character of the relativistic motion has also 
been exposed by the Lattice calculations |^, since they have obtained larger values 


2 



of fs/fo than Md/Mb from the nonrelativistic scaling law. This sitnation can 
be nnderstood clearly throngh onr relativistic calcnlation. 

The potential model has been snccessfnl for and T families with the nonrel¬ 
ativistic hamiltonian, since their heavy qnarks can be treated nonrelativistically. 
However, for D or B meson it has been difhcnlt to apply the potential model 
becanse of the relativistic motion of the light qnark in D or H meson. In onr cal¬ 
cnlation we work with the pnrely relativistic hamiltonian, and adopt the variational 
method 0. We take the Ganssian and hydrogen-type wave fnnctions separately 
as trial wave fnnctions [0, and obtain the gronnd state energy and wave fnnction 
by minimizing the expectation valne of the relativistic hamiltonian. By nsing the 
wave fnnction at origin '0(0), we can obtain the valne of the decay constant fp 
from the Van Royen-Weisskopf formnla. Then we can obtain the ratios of the 
decay constants. The reason why we choose the Ganssian and hydrogen-type trial 
wave fnnctions is that the former one is appropriate to the long range conhning 
potential, and the latter one to the short range asymptotically free potential. 

The heavy-light psendoscalar meson is composed of one heavy qnark with mass 
ttiq and one light qnark with rriq, and its relativistic hamiltonian is given by 

H = -\/p^ + + V^p 2 + + V{r), (1) 

where r and p are the relative coordinate and its conjngate momentnm. The hamil¬ 
tonian in m represents the energy of the meson in the center of mass coordinate, 
since in that reference frame the momenta of both the heavy and light qnarks have 
the same magnitnde as that of the conjngate momentnm of the relative coordinate. 

We apply the variational method to the hamiltonian (J^) with the Ganssian and 
hydrogen-type trial wave fnnctions. The Ganssian wave fnnction is given by 

V.(r) = (2) 

V^r 

where /i is the variational parameter. The Fonrier transform of '0(r) gives the 
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momentum space wave function x(p)) which is also Gaussian, 


x(p) = 




(a/^P)3/2 

The ground state is given by minimizing the expectation value of if in ([^), 
{H) = = 0 at fx = fl, 

dfi 


( 3 ) 


( 4 ) 


and then ft = represents the inverse size of the meson, and E = E{fl) its mass 
Mp p. For the value of the light quark mass rrig in (|l|), we use the current quark 
mass given by Dominguez and Rafael §]: md = 9.9 MeV and rris = 199 MeV. 

We perform the same calculation as the above for the hydrogen-type wave 
function 


V'(r) = 


^-r/ao 


( 5 ) 


where oq is the variational parameter which represents the size of the meson. The 
momentum space wave function conjugate to the '^(r) in is given by 

2^2 


x(p) = 


3/2 

% 


n (agp 2 _^xj 2 ' 


( 6 ) 


For V{r) in (|I]), we consider the following six potentials of the potential model, 
which we also display in Fig. 1. We note in Fig. 1 the tendency that the potential 
which has higher values of potential energy in the short range has lower values in 
the long range, and vice versa. 

(A) Coulomb and linear potential of Eichten et al. Q: 

V{r) = -^ + Kr, (7) 

r 

with etc = 0.52, 77 = l/(2.34)^ GeV^, me = 1.84 GeV, m^ = 5.18 GeV. 


(B) Coulomb and linear potential of Hagiwara et al. 




rv 

V{r) = -^ + Kr, 

r 


( 8 ) 
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with ac = 0.47, K = 0.19 GeV^, rric = 1.32 GeV, mb = 4.75 GeV. 


(G) Power law potential of Martin 
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V{r) = -8.064 GeV + (6.898 GeV) (r ■ 1 GeV)°-\ 


( 9 ) 


with me = 1.8 GeV, mb = 5.174 GeV. 


(D) Power law potential of Rosner et al. \\2 


V{r) = -0.772 GeV + 0.801 ( (r ■ 1 GeV)" - 1 ) / a. 


( 10 ) 


with ot = —0.12, me = 1.56 GeV, mb = 4.96 GeV. 


(E) Logarithmic potential of Quigg and Rosner [|T^ : 


V{r) = -0.6635 GeV + (0.733 GeV) log(r ■ 1 GeV), 


( 11 ) 


with me = 1.5 GeV, mb = 4.906 GeV. 


(F) Richardson potential [|T^ 


with Uf = 3, A = 0.398 GeV, me = 1.491 GeV, 


dq e 

q [ln(g2 — 1)]^ + vr^’ 

mb = 4.884 GeV. 


( 12 ) 


The results of the variational calculations with the Gaussian wave function are 
organized in Table 1, and those with the hydrogen-type wave function in Table 2. 
We see in Table 1 and 2 that the larger energy {E) state has the smaller size of 
meson (l//i or oq). In order to check whether the Gaussian wave function in (^ is a 
really good wave function, we enlarged the trial wave function by adding the second 
excited harmonic oscillator eigenfunction which is an even function, since the hrst 
excited one which is an odd function can not be included in the ground state wave 
function of the relativistic hamiltonian which commutes with the parity operator. 
For this enlarged trial wave function we obtained the result that the Gaussian 
part contributes much dominantly, therefore it conhrms that the Gaussian wave 
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function is a very good trial wave function in the variational calculation of the 
relativistic hamiltonian in with the harmonic type wave function. We also 
checked the wave function in by adding the hrst excited hydrogen-type wave 
function, and conhrmed that the wave function in (|^) is also much dominant for 
the ground state wave function of (|I|). 

The decay constant fp of the pseudoscalar meson P is dehned by the matrix 
element (0|A^|P(g)): 

{0\A,\P{q)) = tqjp. (13) 


By considering the low energy limit of the heavy meson annihilation, we have the 
relation between fp and the ground state wave function at origin 'fp{0) from the 
Van Royen-Weisskopf formula with the color factor [^]: 


/l = ^W.p(o)r. 


(14) 


where Mp is the heavy meson mass. Using this formula, from and we have 

12 /P^(P)n3/2 


/ 22 

fp = J )_ ^ ) ' for Gaussian wave function. 


V Mp V 


TT ^ 

1 \ 3/2 


(15) 


for hydrogen — type wave function. (16) 


"7r^Pao(P)' 

Then the ratio of the decay constants of pseudoscalar mesons A and B is given by 

i 


[a 

Ib 


IM, 


B 


Ma 


X 


( ^ Gaussian wave function, (17) 


= y X ^ ^ ^ for hydrogen — type wave function. (18) 

By using the values in Table 1 and 2 for M (P), 'p^ (fl), and Oq in ([T7|) and 
), we obtain the ratios of the decay constants for the potential models (A)- 
(F), which we present in Table 3 and 4. We see in Table 3 and 4 that fsAfos 
is enhanced, compared with the nonrelativistic scaling law \JMp,^ /Mp^ whose ex¬ 
perimental value is 0.605, by the factor of 1.324 which is induced by the factor 
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of |V’-Bs(0)/'^i5s(0)l- For the estimation of this enhancement factor we used the 
average of the fourth columns of Table 3 and 4. Sometimes this factor has been 
approximated to be 1 and /bs/Ids — \J has been used, by treating it in 
analogy with the nonrelativistic case |^. However our calculation shows that this 
factor is indeed important and different from 1 signihcantly. The same situation 
happens for the ratio /b^//u)^. Compared with Mb^/Mb^ whose experimental 
value is 0.595, fBa/foa is enhanced by the factor of 1.335, which is estimated from 
the average of the fifth columns of Table 3 and 4. This situation is in agreement 
with the results of the Lattice calculations p, as we see in Table 5, where we 
organized the results of Lattice calculations. 

As an application of the obtained ratios of decay constants, let us consider 
the determination of the values of fBsi ^-nd Jb^ from the experimental value 
of fos- The cleanest way to obtain the value of the decay constant from the 
experimental results is through the purely leptonic decays of the meson, which 
are understood theoretically to occur via an annihilation of the two valence quarks. 
The decay rate of the meson is given by the formula 

1 ■ -2 


r(Z),+ ^ ^ Gl fl, m? Mo. (l 


mf 


IV 


2 

CS I ? 


(19) 


where /d„ is the meson decay constant, Mbi^ is the Dg mass, mi is the mass of the 
hnal-state lepton, Gp is the Fermi coupling constant, and Vcs is the CKM matrix 
element. The WA75 and CLEO collaborations took the data for the branching 
ratio and the Review of Particle Properties presents 

B{Df —>■ = (5.9 ± 2.2) X 10“^. By using this branching ratio, the life time of 

Dg meson r = (0.467 ± 0.017) x lO-^^ ^ ^ggg ^ ^ g y ^eV, = 105.66 

MeV, Gp = 1.1664 x 10“® GeV“^, and \Vcs\ = 1.01 ± 0.18 [^|, we obtain the 
following value of from (0): 


/d, = 265 ± 68 MeV. 


( 20 ) 


The uncertainty of the value of is due to those of the experimental values of 
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B{Df —>■ l^v) and \Vcs\^ therefore if their values are improved experimentally, we 
can obtain very accurately. When we combine the value in (pOf ) and the 
ratios given by the average of those in Table 3 and 4, we get the following values of 
/b,, and fs/. fn, = 253±65 MeV, fs, = 212±54 MeV, fs, = 201±51 MeV. 

In summary, we calculated the various ratios of the B and D meson decay 
constants by applying the variational method to the relativistic hamiltonian. We 
took the Gaussian and hydrogen-type trial wave functions separately, and used six 
different potentials for the potential energy. We obtained the results which are 
similar for different trial wave functions and different potentials. This fact implies 
that our method is reliable. The obtained results for the ratio fs/fo show that 
the nonrelativistic scaling law fs/fo — \JMd/Mb should be implemented by the 
relativistic consideration. Its enhancement factor we obtained is about 1.33, which 
is induced by |V’n(O)/'0D(O)|. This result is in a pretty good agreement with the 
recent Lattice calculations. Our results for fBs !and fos/foa are both about 
1.05. We also obtained the value of the double ratio {fBs/fBa)/{fDs/fDd) as 1.008, 
whereas Grinstein obtained 0.967 with the heavy quark effective theory, and Oakes 
1.004 with the chiral symmetry breaking. 
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Model 

KBs) 

E{B,) 

KBd) 

E{Ba) 

KBs) 

E{Ds) 

KDd) 

E{Dd) 

A ( Eich. ) 

0.565 

5.933 

0.544 

5.896 

0.495 

2.661 

0.478 

2.620 

B ( Hagi. ) 

0.550 

5.553 

0.530 

5.516 

0.466 

2.217 

0.451 

2.174 

C (Power 1) 

0.587 

5.154 

0.567 

5.119 

0.510 

1.855 

0.493 

1.815 

D (Power 2) 

0.600 

5.353 

0.577 

5.318 

0.495 

2.042 

0.475 

2.000 

E ( Log. ) 

0.586 

5.396 

0.565 

5.360 

0.490 

2.080 

0.472 

2.038 

F ( Rich. ) 

0.584 

5.404 

0.564 

5.368 

0.495 

2.103 

0.479 

2.062 

(Average) 

0.579 

5.466 

0.558 

5.430 

0.492 

2.160 

0.475 

2.118 


Table 1: The values of the variational parameter fi which minimize (H), and the 
corresponding values of the minimum energy for the Gaussian wave function 


Model 

o-oiBg) 

E{Bg) 

o-oiBd) 

E{Ba) 

o-o{Ds) 

E{D,) 

Oo(-Dd) 

E{Da) 

A ( Eich. ) 

1.260 

5.884 

1.310 

5.845 

1.549 

2.140 

1.605 

2.094 

B ( Hagi. ) 

1.309 

5.512 

1.360 

5.472 

1.570 

2.189 

1.625 

2.143 

C (Power 1) 

1.305 

5.118 

1.353 

5.078 

1.504 

1.823 

1.561 

1.779 

D (Power 2) 

1.264 

5.303 

1.317 

5.264 

1.532 

1.997 

1.602 

1.951 

E ( Log. ) 

1.308 

5.200 

1.360 

5.160 

1.593 

1.876 

1.660 

1.829 

F ( Rich. ) 

1.260 

5.359 

1.307 

5.321 

1.503 

2.068 

1.558 

2.024 

(Average) 

1.284 

5.396 

1.335 

5.357 

1.542 

2.016 

1.602 

1.970 


Table 2: The values of the variational parameter oq which minimize (H), and the 
corresponding values of the minimum energy for the hydrogen-type wave function. 
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Model 

fsj fsa 

JdJ fDa 

UBs/ fB^) 

/d,//b. 

fod/fsd 

A ( Rich. ) 

1.055 

1.046 


1.225 

1.236 

B ( Hagi. ) 

1.054 

1.040 


1.234 

1.250 

C (Power 1) 

1.050 

1.041 


1.350 

1.362 

D (Power 2) 

1.057 

1.053 


1.213 

1.218 

E ( Log. ) 

1.053 

1.047 


1.232 

1.238 

F ( Rich. ) 

1.050 

1.040 


1.251 

1.263 

(Average) 

1.053±.003 

1.045±.005 

1.008±.003 

1.251±.046 

1.261±.047 


Table 3: Ratios of the decay constants obtained for the Gaussian wave function. 


Model 

JbsI fSd 

fOs/fOd 

{Ibs/Ib^} 

(fOs/fOd) 

Jds/ fss 

fod/f Bd 

A ( Eich. ) 

1.057 

1.043 


1.217 

1.232 

B ( Hagi. ) 

1.055 

1.042 


1.208 

1.223 

C (Power 1) 

1.052 

1.044 


1.354 

1.363 

D (Power 2) 

1.060 

1.057 


1.221 

1.224 

E ( Log. ) 

1.056 

1.050 


1.239 

1.246 

F ( Rich. ) 

1.053 

1.044 


1.236 

1.246 

(Average) 

1.056±.003 

1.047±.005 

1.008±.004 

1.246±.050 

1.256±.049 


Table 4: Ratios of the decay constants obtained for the hydrogen-type wave func¬ 
tion. 
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Group 

IBs / fBa 

fOs/f Da 

(fOs/fOa) 

IDs / IBs 

fOd/f Bd 

ELC PI 
UKQCD |16 
BLS |17| 
MILC 18 

1.08±.06 

1 99 +.04 
-.03 

1.11±.02±.05 

1.13(2)(9)(4) 

1.08±.02 

1.18±.02 

1.11±.02± .05 
1.09(1)(4)(4) 

l.OOi.06 

1.03 

1.00±.03± .06 
1.04(2)(9)(5) 

1.03±.22 

1 nQ+-04+-42 
-■-•Uy -.03-.06 

1.11±.06± .27 
1.18(3)(17)(13) 

1.02±.21 
lie +.05+.46 
-‘-•J-U _,05-.l4 

1.11±.08±.30 

1.22(5)(17)(19) 


Table 5: The results of the ratios from Lattice calculations: the third column was 
estimated from the values of the hrst and second columns, and the fourth and hfth 
columns were estimated from the Lattice calculation results of /d^, and 

fBd- 




























horizontal axis is in the unit of GeV ^ (1 GeV ^ = 0.197 fm), and the potential 
energy of the vertical axis is in the unit of GeV. 
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